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We study the quantum phase transition in an atomic Bose gas near a Feshbach resonance in terms 
of the renormalization group. This quantum phase transition is characterized by an Ising order 
parameter. We show that in the low temperature regime where the quantum fluctuations dominate 
the low-energy physics this phase transition is of first order because of the coupling between the Ising 
order parameter and the Goldstone mode existing in the bosonic superfluid. However, when the 
thermal fluctuations become important, the phase transition turns into the second order one, which 
belongs to the three-dimensional Ising universality class. We also calculate the damping rate of the 
collective mode in the phase with only a molecular Bose-Einstein condensate near the second-order 
transition line, which can serve as an experimental signature of the second-order transition. 

PACS numbers: 05.30.Jp, 64.60.Ak, 67.40.Db 



I. INTRODUCTION 

Trapped dilute cold atomic gases are one of the most 
exciting fields in condensed matter physicsiSi An impor- 
tant recent development in this area is the application 
of Feshbach resonances. A Feshbach resonance in the 
scattering amplitude of two atoms occurs when the total 
energy of the atoms is close to the energy of a molecular 
state that is weakly coupled to the atomic continuum. 
Especially, the energy difference between the molecular 
state and the two-atom continuum, known as the detun- 
ing 5, can be experimentally tuned by means of a mag- 
netic field. Therefore, by sweeping the magnetic field 
from positive to negative detuning through the Feshbach 
resonance, it is actually possible to form molecules in the 
atomic gasi^i^ In fact, recently, it has been possible to 
create a Bose-Einstein condensate (BEC) of molecules in 
an atomic Fermi gas with a Feshbach resonancei^ii^ This 
offers the opportunity for observing a Bardeen-Cooper- 
Schrieffer (BCS) transition in a dilute Fermi gas. 

In the present paper, we study an analogous situation 
by varying S in an atomic Bose gas with a Feshbach reso- 
nance. Two recent worka^iifl have shown that by varying 
d there is a true quantum phase transition (QPT) in an 
atomic Bose gas, in contrast to the case of an atomic 
Fermi gas where a smooth BEC-BCS crossover exists as 
one changes S. An argument to understand this has been 
given in Ref. ^3 We recapitulate it in the following to fix 
our notation. The effective Lagrangian describing a di- 
lute atomic gas with a Feshbach resonance can be written 
as 



— a ci Q' Q' ' — " ' 



(1) 



operators of atoms and molecules, respectively, m and M 
are the mass of the atom and that of the molecule, respec- 
tively. In the path integral formula, \E'm is an ordinary 
number. On the other hand, is an ordinary number 
for the atomic Bose gas, while it is a Grassmann number 
for the atomic Fermi gas. (For the atomic Fermi gas, 
contains the indices for hyperfine spins and should be un- 
derstood as a spinor.) The Lagrangian L [Eq. |^] has a 
U(l) symmetry. That is, it is invariant against the U(l) 
transformation 
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2 ati^^ 2 
with /i = 2/i — (5. Here ^'a and arc the annihilation 



For an atomic Bose gas, due to the go term, a nonzero 
value of {'^a) (atomic BEC) must lead to a nonzero value 
of (^'m) (molecular BEC). However, the reverse is not 
true. That is, it is possible for the gas to contain only 
a molecular BEC. In the case with both atomic BEC 
and molecular BEC, the U(l) symmetry of L is com- 
pletely broken. But for the case with only molecular 
BEC there is a residual Z2 symmetry. That is, the effec- 
tive Lagrangian in this case is invariant against the Z2 
transformation: — ^'a and — > ^m- The above 

analysis indicates that at low temperature the atomic 
Bose gas with a Feshbach resonance has two thermody- 
namically distinct phases: the "atomic superfluid" (ASF) 
phase with both atomic BEC and molecular BEC, and 
the "molecular superfluid" (MSF) phase with molecular 
BEC only. For an atomic Fermi gas, due to the 5^0 term, 
a nonzero value of {^a'^a) must be accompanied with a 
nonvanishing value of {^m) and vice versa. Therefore, 
the BCS region has the same symmetry as the BEC re- 
gion, and only a crossover occurs. 

Based on the above observation, the MSF phase is dis- 
tinct from the ASF phase by a Z2 symmetry. Thus, the 
phase transition between the two phases is characterized 
by an Ising (Z2) order parameter. Naively, one may ex- 
pect that this QPT is of second order and belongs to 
the four-dimensional (4c?) Ising universality class at zero 
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FIG. 1: The schematic phase diagram of the dilute Bose gas 
with single Feshbach resonance near the transition line be- 
tween the ASF and MSF phases below Tc{5), where Tc{5) rep- 
resents the transition temperature between the normal Bose 
gas and the superfluid phase. The dotted line denotes the 
crossover from low temperature to high temperature regimes. 
The solid line is the phase boundary between the ASF and 
MSF phases. The point C is a tricritical point. On the phase 
boundary, the portion below the point C is of first order while 
the one above the point C is of second order. 



temperature and the three-dimensional (3(i) Ising univer- 
sality class at finite temperature. However, there is a gap- 
less excitation in the MSF phase, which is the Goldstone 
mode associated with the molecular BEG. An additional 
gapless excitation may result in severe IR divergences and 
modify the critical behavior. Therefore, a proper the- 
ory describing the QPT should consist of the Ising order 
parameter as well as the Goldstone mode. A renormal- 
ization group (RG) analysis indicates that the coupling 
between the Ising order parameter and the Goldstone 
mode will drive the zero-temperature phase transition 
to become weakly first order— through the Colemann- 
Weinberg mechanism— The question in which we are 
interested here is the nature of this phase transition at 
finite temperature. We study this problem in terms of 
the RG. Our main results are shown in Fig. ^ which are 
valid at the temperature much lower than the amplitude 
fiuctuations in the molecular BEG. At low temperature 
where the quantum fluctuations dominate the low-energy 
physics, the phase transition between the MSF and ASF 
phases is of weakly first order. When the thermal fluctu- 
ations become important, the effects of the coupling be- 
tween the Ising order parameter and the Goldstone mode 
are suppressed and the transition becomes a second-order 
one belonging to the id Ising universality class. Thus, a 
tricritical point must exist on the phase boundary to sep- 
arate these two kinds of phase transition. 

The rest of the paper is organized as follows: We write 
down the effective Lagrangian describing the QPT by 
symmetry argument and perform one-loop RG analysis in 
Sec. m The solution of the scaling equations is presented 
in Sec. IIIII In Sec. IIVI we calculate the damping rate 



of the collective excitation in the MSF phase near the 
second-order transition line. The last section is devoted 
to our conclusion. 



II. LANDAU THEORY AND 
RENORMALIZATION GROUP ANALYSIS 

A. Effective Lagrangian 

We start with the "normal phase" for this QPT. The 
fundamental fields of the effective theory describing the 
QPT consist of the Ising order parameter 0, which can 
be taken as the imaginary part of '^a, and the phase 
of ^rai d. Following Landau, the corresponding effec- 
tive Lagrangian can be written down through symmetry 
consideration. The symmetries involved here are the Z2 
symmetry, a subgroup of the U(l) symmetry [Eq. 0], 
under which the A and 6 fields transform as 



9-^9 + 27: 



(3) 



and the U(l) symmetry under which 6 transforms as 
9{x) — > 9{x) + a where a is an arbitrary constant. The 
most general local action consistent with the Z2 and U(l) 
symmetries is J = /0 + + lint where 



drdPx ■ 



with uo > 0, 
le = 

and 
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Here to measures the (mean-field) distance from the tran- 
sition point, {to > in the MSF phase while to < in 
the ASF phase.) ps is the "bare" superfluid density and 
Co and vq are "bare" velocities for the order parameter 
(j) and the Goldstone boson 9, respectively. In Eq. ((BJ, 
only the most relevant (near the transition point) cou- 
pling between the Goldstone mode and the order param- 
eter is kept. The factor i in Eq. © is dictated by the 
charge conjugation symmetry of the original Lagrangian 
L [Eq. (UJ]. It requires that L ^ when *1 and 

'^m ^\n^ which leads to / ^ when (f) —0 and 
9 —9. The natural cutoff of this action is provided by 
the gap of amplitude fluctuations in the molecular BEG. 
A similar action has appeared in other context— 

The action Q — ^ can also be derived from the 
Lagrangian L [Eq. J^l] by considering the fiuctuations 
around its mean-field solution in the MSF phase. By in- 
tegrating out the sectors with finite gaps at the transition 
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point, i.e. the real part of ^'a and the amphtude of 

in terms of the perturbation theory in ui, U2, M3, and 

go/y/n, to the leading order, one may obtain 
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near the transition point. (Here we take the real part 
of '^a as the gapped sector at the transition point. This 
amounts to assuming go > 0.) In Eq. ((TJ, n is the 
total particle density and i5o denotes the magnetic field 
detuning at which the zero-temperature quantum phase 
transition occurs at the mean-field level, which is given 
by 5o « (u3 - ^) n - 23/2^0^1/2. (This value of So has 
already been given by Ref. From Eq. (7J), it is ob- 
vious that the value of uo may become negative in some 
parameter regime. For uq < 0, one should include the 
terms with higher powers of (f> like (j>^ in to make it 
stable, and the QPT between the MSF and ASF phases 
will be of first order ^ even in the absence of the coupling 
to the Goldstone mode 6. In this case, it is of little in- 
terest because the correlation length is always finite and 
thus there are no universal behaviors for physical quan- 
tities at both zero and finite temperatures. Therefore, 
in the following, we shall concentrated on the parameter 
regime which gives rise to a positive value of uq. From 
Eq. 0, this corresponds to (M3 -|- 1/2/71 go)^ < UiU2- 

AtT ~ 0, with Uq > describes the M Ising univer- 
sality class. It exhibits a Gaussian behavior with loga- 
rithmic corrections due to the presence of the marginally 
irrelevant coupling uq. Therefore, we may consider the 
RG transformation 



t/.s , X ^ x/s 
sch , s0 , 



(8) 



while Co, Uq, Ps, and M remain invariant. Here s > 1 
is the rescaling factor. Then, the action + Ig + lint 
with mq = = Aq is invariant under the RG transfor- 
mation [Eq. JHJ]. It is straightforward to see that the 
coupling constants uo and Aq are marginal at the tree 
level with respect to the Gaussian fixed point under the 
RG transformation [Eq. ©]. Within the weak-coupling 
region, one-loop RG equations are needed to determine 
their roles on the low-energy physics. 



B. One-loop RG equations 

To compute the one- loop RG equations, we employ 
the momentum-shell RG. Before integrating out the fast 
modes, we make a change on the variables: 



X ^ A ^x 



Ac, 



1/2, 



A Cq T , w = vo/co 
1/2 



r ^ A ^to , u = cqWo , P 



Acq 
T 



where A is an UV cutoff in momenta and T is the tem- 
perature. Then, v, r, u, A, and /3 become dimension- 
less parameters and our working action is written as 



/ = /q J(Px C where 



1 



dr0] +(Ve)2 
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(9) 



By integrating out the fast modes, i.e. those modes 
with momenta within the momentum shell e~' < |fc| < 1 
where e~' is the scaling factor, we obtain the one-loop 
RG equations: 



dt. 



Km 



dl 

^ = -4(i+V>/2 /2(x/TT7?A0 

dri UiU 
dl '^4(l + rO 



— ^ gi(VlTnltuvi/ti) 
ti 



(10) 

(11) 

(12) 



dKi 
dl 



4 (l + r;)3/2 

2Kf 



VI + n(l + ri- vf) 



.g2(\/l + n/ti,vi/ti) 
hivi/ti) , (13) 



dl 



-24K 



2'iKf 
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Vl 

(1 + n -1)2)2 

6KIUI 



53(Vl + ri/ti,vi/ti) 

.94 (Vl + n/tuvi/ti) 
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where t = \/l3,K= (Au)V(27r2), U = uI{2t{^) and 
~ a-coth( 



/2(a;)=coth(|) +|sinh 
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Here the rescaled quantities have exphcit I dependence 
indicated, e.g. t/, while quantities without I dependence 
(e.g. T) refer to physical quantities. 

Since K > 0, we always have dvi/dl < from Eq. 
and thus vi will flow to zero. Therefore, we may set 

= in Eqs. I^^ll — yielding 



dri 
If 


- 2ri 

1 + n 


, iUi + 
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tl 



The solution of Eq. (|10|1 is simply given by 
tl = teJ . 



(15) 
(16) 
(17) 

(18) 



In the next section, we shall solve Eqs. p5|l — (|17|) ap- 
proximately by neglecting the possible logarithmic cor- 
rections. 



III. SOLUTION OF SCALING EQUATIONS 

Scaling stops when |n| ~ 1. One must distinguish two 
regimes: ti <^ 1 and 3> 1. The former corresponds to 
the quantum {T = 0) regime where the quantum fluctu- 
ations dominate the low-energy physics. Otherwise, it is 
the continuum high T regime where the thermal fluctu- 
ations become important. 



A. Quantum regime 

In the quantum regime, the nonlinear dependence of 
the RG functions (those terms appearing at the R.H.S. of 
Eqs. (|15|) — (|17|l ) on r; will not be important. Thus, we 
shall consider the limiting equations, valid for r, K^U <^ 
1 
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FIG. 2: Typical solutions of Eqs. with the initial 

conditions (a) r = 0.002, K = 0.0003, and U = 0.005, (b) 
r = 0.002, K = 0.003, and U = 0.005, and (c) r = 0.002, 
K = 0.03, and U = 0.005. The x-axis is r; and the y — axis 
is Ki/K or Ui/U. 



(20) 



^ = 24i^Kt^; + + 8A'0'/2(^^) • (21) 

To obtain the condition on T for the occurrence of the 
quantum regime, one may set t = in Eqs. p9|l — 121|l . 
yielding 



dl 6 12 



dl 
dyi 



1 



dXi 



(3^i-4)(3;,-l) . 



(22) 
(23) 
(24) 



Here Xi = Ui + l2Ki and yi = Ui/Xi. The solution of 
Eqs. 123) and is written as 



X, 



Ui 



Xq 



1 + 3xol/8 



doXl^^ - 4 



doXi 



1/3 



(25) 



where xq = U + 12K and the value of do is chosen such 

1 /3 

that Ui — U when Xi — x^. (In fact, doa^o = 4 + 
U/{AK). Therefore, do > 1 for K, U < 1.) 

By solving Eqs. I|22() — H24() numerically, we notice 
that the value of Ki remains almost unchanged when 
ri < 1, as shown in Fig. |21 However, the behavior of Ui 
depends on the sign of the function G{r,U, K), roughly 
given by 



G(r, t/, K] 



16 

3x0 



U 
T6K 



ln(l/r). (26) 
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(Equation H26() is determined from Eq. H25|l by requiring 

1 /3 

that Ui* = or doX^J =4 where n. = 1. The con- 
dition r/. = 1 roughly gives rise to 21* w ln(l/r).) (i) 
As G{r,U,K) > 0, [/; > and its magnitude does not 
change much for n < 1, as shown in Fig. [21(a) and (b). 
(ii) As G{r, U, K) < 0, the vahic of Ui becomes negative 
rapidly before r; = 1 as shown in Fig. [21 (c). 

For given K and U, we may increase the value of r 
starting with r = 0. In the beginning, G{r, K,U) < 
and U becomes negative when scaling stops. Thus, a 
(j>^ term must be included in the effective action. A 
mean- field theory in this case gives rise to {(p) ^ 0. 
This corresponds to the ASF phase. On the other hand, 
for sufficiently large r such that G{r,U,K) > 0, t7 is 
still positive when scaling stops, and thus a mean-field 
theory gives rise to {(f>) = 0. This corresponds to the 
MSF phase. Consequently, there are two phases at zero 
temperature separated by a transition point, denoted by 
fco- The value of Tco can be estimated by requiring 
that G{rco,U, K) = 0. A negative value of U is usu- 
ally interpreted as a fluctuation-induced first-order phase 
transitioniii Therefore, the QPT at T = is of first or- 
der. In the following, we shall focus on case (i), i.e. the 
Z2 symmetric (MSF) phase. 

In case (i), because the values of Ki and Ui do not 
change significantly when scaling stops, we may set Ki ~ 
JC > and Ui =U > in Eq. ((221) for simplicity, where 
IC/K w 1 and U/U < 1. Within this approximation, the 
logarithmic corrections are neglected. (The approxima- 
tion we use amounts to neglecting the / dependence of 
Xi. Inclusion of it will result in logarithmic corrections.) 
Setting n = 1, solving for substituting the result into 
Eq. H18|l . and demanding <C 1, one may obtain the 
condition for the occurrence of the quantum regime in 
the MSF phase: 



t < 



1 + Ai 



1/2 



(27) 



where Ai = JC+U/%. 

In terms of Eq. Q, the crossover line between the 
quantum and continuum high T regimes in the MSF 
phase, corresponding to r > Tco, can be expressed by 



h + 



coA2 
27r2 



U2\ 



Uq 

8 



(28) 



as shown in Fig. ^ Here the values of cq, to, Ao, and uq 
near the transition point is given by Eq. ((TJ. In addition, 
the momentum UV cutoff A can be estimated by A « 
■min{2^^^y/mgo n^/^, ^uiMj^ n^/^}, where the former 
and the latter arise from the gap of the real part of '^a. and 
that of the amplitude fluctuations of <f>„i, respectively. To 
obtain Eq. (|^, we have neglected 1 -f Ai in Eq. lf77|) 
because Ai ^ 1. Further, K, and U are replaced by K 
and ?7, respectively. 

Finally, we determine the correlation length in the 
quantum regime. To do so, we may set Ki = K, and 



Ui =14 in Eq. H19|l and performing low-temperature ex- 
pansion, yielding 



^ = 2ri+ 2Ai - AJCte^ . 
dl 



(29) 



The solution of Eq. H29|l is written as 

n = e^' [r + Ai (1 - e-2') - (l - e"')] • (30) 

The correlation length is given by ~ where = 1. 
From Eq. (|30|) . one may obtain 



^^{r + A^-A^t)-^'^ 



(31) 



with A2 = 4/C. Therefore, ^ is always finite for given 
t satisfying the inequality H27|) . as expected for a first- 
order phase transition. To sum up, the phase transition 
between the MSF and ASF phases in the quantum (low 
temperature) regime is of weakly first-order. 



B. Continuum high T regime 

In this regime, it is convenient to divide the scaling 
into two steps: ti -^1 and ^ 1. This introduces mul- 
tiplicative errors of order unit coming from the imprecise 
treatment of the crossover regime ti ^ \. In the first 
step, we integrate over I from to ^, such that = 1. 
Next we consider the case t; 3> 1. In this case, it is more 
convenient to define the coupling V/ = tiJJi. We would 
like to show that the RG equations for r/ and Vi in the 
continuum high T regime are identical to those for the 
3d Ising model. 

We start with Eqs. and and take > 1. 

Using /i(0) = 2 and f2{x ^ 0) w 4/a;, we obtain 



±i = 2n+ ^' , 

dl ' 2(l-f rO ' 

dl ' 2{\ + riY ' 



(32) 
(33) 



Equations l|32() and (|33|l are identical to the one-loop RG 
equations for the 3c? Ising model. Therefore, we conclude 
that the transition between the MSF and ASF phases 
in the continuum high T regime is of second order and 
belongs to the ?>d Ising universality class. Near the tran- 
sition point, the solution of Eqs. (|32|) and 133|) is given 

by 



(34) 



where vj, w 0.63 is the correlation length exponent for the 
3(i Ising model and f = rj, which are determined by the 
RG equations in the quantum regime. Here r* and V* are 
the fixed points of Eqs. H32|) and H33|) . which within the e 
expansion corresponds to the Wilson-Fisher fixed point, 
where e = 4 — c? and d is the spatial dimensions. The ex- 
act values of r* and V* arc nonuniversal, and the deter- 
mination of them is beyond the present field-theoretical 
approach. 
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To determine r, we shall focus on the MSF phase. In 
this case, one may neglect the possible logarithmic cor- 
rections and set Ki ~ K- > and Ui ~ U > where 
IC/K « 1 and U/U < 1. (Note that both IC and U may 
not be equal to IC and U.) Then, from Eq. (|19|l . we 
obtain 



Bi 



t2 



B2 
t 



c 



(35) 



where Bi = K.+U/?>, B2 = 4/C, and 



4Si / dy y-^'uBiy-^) , 
/o 



with nsix) = (e^ - 1)"^ 
Scaling stops when |n 

length is given by ^ 
gives rise to 



r 



-j/aln 



Si 



The correlation 
Equation 



i2 



B2 
t 



B. 



where i?3 = C — r* is a nonuniversal constant. As a 
result, we get 



calculate the damping rate of the collective mode in the 
MSF phase near the second-order transition line, arising 
from the coupling to the critical fluctuations. Since the 
identification of the Ising transition is established within 
the e expansion, the following RG-improvcd perturbative 
calculation is supposed to be understood within such a 
context. 

The full propagator of the 9 field can be written as 
D~^{iuJn,k) ~ DQ^{iu!mk) + T,g{iujmk) where aj„ = 
2n7r//3, Se is the self-energy of the 9 field, and 



Do{k) = 



1 



is the free propagator of the 9 field. The spectrum is 
determined by 



D ^ {iuj„ 



iO+,k) = , 



which gives rise to w = — ijk where 7^ denotes the 
damping rate and Cfc is the dispersion relation of the col- 
lection mode. In the long wavelength limit, = vk 
where k = \k,\ and v is the dimcnsionlcss rcnormalized 
velocity measured in c. By solving the equation, one 
may obtain 



\r-rc{t)\ 



(36) 



where rdt) = —Bi + — B^t^ denotes the transition 
point between the MSF and ASF phases for given tem- 
perature. Equation H36|l is valid up to logarithmic cor- 
rections. 

Finally, we consider the RG equation for Ki in the 
continuum high T regime, which can be written as 



Ui 



dRi 



R 



I 



Ri 



dUi 3 3 
where R = K/U . Solving Eq. 1)3 7|) gives rise to 

V* , 



(37) 



(38) 



near the transition line. When scaling stops, wc may 
perform the perturbation theory in Ki and Vi as long 
as = V"*/(4^i) < f. Within the spirit of the e 

expansion, it is indeed the case because V* = 0(e). 



IV. DAMPING RATE OF COLLECTIVE 
EXCITATIONS IN THE MSF PHASE 

The damping in the collective excitations of conden- 
sates will become much more severe near the second- 
order transition line due to the coupling to the critical 
fluctuations. Therefore, a measurement of the damping 
in the collective modes may serve as a signature of the 
second-order Ising transition between the MSF and ASF 
phases. In this section, we shall use the RG equations to 



7fc = — lmT:e{(^ + «0+, k) 



(39) 



with the understanding that lo is replace by ek- Near the 
second-order transition line, = F(k, K, U, t, A) where 

A = r — r* and we have 



F{k, K, U, t, A) = e-'F{e'k, Ki, Ui,e% A,) , 

because the scaling dimension of 7fc is [7fc] — 1. Scaling 
stops at A;* = 1, and wc obtain 



Ik 



(40) 



where ^{xi, X2, X3, X4) — F(xi,X2,X3,X4,l) is the scal- 
ing function. 

To the one-loop order, the scaling function $ can be 
written as 

$(^fc,A>,C/,.,ei) 

Ki' f^3^ S[e^k - E{p) + E{^k + p)] 



E{p)E{^k+p) 



X < UB 



E{p) 



Ki*e^k 



16 



X < UB 



d^p 



- riB 


\E{^k+p^ 


} 


[ \ 




S[eik 


-E{p)-E{ik+p)] 



E{p) 



- riB 



E{p)E{Ck+p 
E{^k+p 



Ki-€ik [,3 S[e^k + E{p) + E{^k+p)] 



16 



X < UB 



E{p) 



riB 



E{p)E{^k+p 
E{^k+p 



(41) 
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Here E{p) = yjip^ + 1 and nsix) = (e^ — 1)"-^. For small 
k, we find that the damping rate to the lowest order in k 
is given hyi^ 



Ik 



nV*v^ / k 
16 W 



1 



exp 



-1 



(42) 



where from Eq. ^ ^ [\r - rc{t)\/t'^] From Eq. 
()42|l . we see that away from the critical region, the damp- 
ing of the collective excitations arising from the coupling 
to the critical fluctuations is exponentially small. On the 
other hand, diverges near the transition line, and thus 
within the critical region, jk can be approximated as 



7fe 



16 



v'k + 0[ ^ 



(43) 



Therefore, the collective mode is heavily damped within 
the critical region such that no well-defined collective ex- 
citations exist there on account oijk/ek — 0(1). 

Due to the coupling to the amplitude fluctuations, the 
collective mode in the BEG already acquires a damping 
rate which exhibits dramatic temperature dependence as 
given by Ref. Il6i For example, in the low temperature 
limit, it is of the form: 7^ cx T'^ki^ Compared with 
that result, the damping rate due to the coupling to the 
critical fluctuations in the critical region is indeed much 
stronger than the one arising from the coupling to the 
amplitude fluctuations. However, a careful measurement 
of the damping rate near the second-order transition line 
must be conducted to extract the 3d Ising correlation 
length exponent 1^3. Because within the critical region 
the damping is insensitive to the variation of temperature 
as shown by Eq. H43|l . 



V. CONCLUSION 

We study the QPT between the MSF and ASF phases 
in terms of RG. Our calculations suggest that in some 
parameter region this transition is of first order as long 
as the temperature is below the transition temperature 



from the normal Bose gas to the supcrfluid phase. How- 
ever, in the other parameter region, the physics of the 
QPT is more interesting. It is of weakly first-order in 
the low temperature regime where the inequality given 
by Eq. H27() is satisfied. As the temperature is raised 
to the continuum high T regime where the inequality 
given by Eq. H27|) is reversed, the transition becomes 
a second-order one, belonging to the 3d Ising universal- 
ity class. This second-order transition in the continuum 
high T regime is guaranteed by two facts: (i) Those terms 
containing the coupling constant K in the one-loop RG 
equations exactly cancel each other in the high temper- 
ature limit, (ii) The coupling constant K becomes an 
irrelevant operator around the 3d Ising (Wilson-Fisher) 
fixed point in the continuum high T regime and flows to 
zero. A schematic phase diagram is shown in Fig. ^ and 
the crossover line between the quantum and continuum 
high T regimes is expressed by the parameters which are 
directly related to experimentally measurable quantities, 
as given by Eq. We must emphasized that all our 

results are valid only at the temperature much lower than 
the gap of amplitude fluctuations in the molecular BEG. 

We also calculate the damping in the collective excita- 
tion in the MSF phase near the second-order transition 
line. Because of the coupling to the critical fluctuations, 
the damping will be strongly enhanced, which may serve 
as an experimental signature of the second-order tran- 
sition. Our calculation shows that the damping rate of 
the collective mode in the MSF phase is insensitive to 
the variation of temperature near the second-order tran- 
sition line, a behavior very different from the one in the 
ordinary BEG where the damping is much weaker and 
exhibits dramatic temperature dependence. 
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